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We have analyzed the beam spin asymmetry and the longitudinally polarized target spin asymmetry of 
the Deep Virtual Compton Scattering process, recently measured by the Jefferson Lab CLAS collabora- 
tion. Our aim is to extract information about the Generalized Parton Distributions of the proton. By fitting 
these data, in a largely model-independent procedure, we are able to extract numerical values for the two 
Compton Form Factors Hj m and Hj m with uncertainties, in average, of the order of 30%. 
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The study of Generalized Parton Distributions (GPDs) 
is currently one of the most intense fields of research in 
hadronic physics, theory-wise as well as experiment-wise. 
GPDs give access in an unprecedented way to part of the 
complex composite structure of the nucleon (or more gen- 
erally of hadrons), which, until now, is not fully calcula- 
ble from first principles of Quantum Chromo-Dynamics 
(QCD). For instance, nucleon GPDs encode, in a frame 
where the nucleon has a quasi-infinite momentum in a cer- 
tain direction (the so-called "infinite momentum frame"), 
the longitudinal momentum distributions of the quarks and 
gluons in the nucleon, their transverse spatial distribution 
and, overall, the correlation between these two latter dis- 
tributions, which is brand new information. As a con- 
sequence of this longitudinal momentum-transverse space 
correlations, there is the possibility to access the contri- 
bution of quarks to the orbital momentum of the nucleon. 
This is of great interest for the notorious "spin puzzle" 
of the nucleon, a long-standing issue in nucleon structure 
studies. We refer the reader to Refs. OMi for the origi- 
nal theoretical articles and recent comprehensive reviews 
on GPDs and for details on the theoretical formalism. 

Nucleon GPDs are the structure functions which are ac- 
cessed, through the factorisation property of QCD, in the 
hard exclusive electroproduction of a meson or a photon 
off the nucleon. If we focus on quark GPDs, the golden 
channel to access them is the Deep Virtual Compton Scat- 
tering (DVCS) process, due to the purely electromagnetic 
nature of the perturbative part of the "handbag" diagram. 
This latter diagram is schematized in Fig. [T] At large 
Q 2 = (e! — e) 2 and small t = (p — p 1 ) 2 , this process 
in which the same quark (or antiquark) absorbs the in- 
coming virtual photon and radiates the final real photon, 
is predicted to be the dominant one. The quantities x + £ 
and x — £ are the longitudinal momentum fractions of the 
initial and final quark (or antiquark) respectively, where 

i = 2=fe and x b = £ (with v = E e i E e ) is the 
standard Deep Inelastic Scattering (DIS) variable. 

For DVCS on the proton, several experimental observ- 
ables measured in different kinematical regimes have been 
published this past decade : cross sections (beam-polarized 
and unpolarized) from the JLab Hall A collaboration 
beam spin asymmetries (BSA) ifiol [Till and longitu- 
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FIG. 1: The handbag diagram for the DVCS process on the proton 
ep — } e'p'j' . There is also a crossed diagram which is not shown 
here. 



dinally polarized target asymmetries (1TSA) [12] from 
the CLAS collaboration and a series of correlated beam- 
charge, beam-spin and transversely polarized target spi n 
asymmetries from the HERMES collaboration Il3l4l5ll . 
The question arises: from this large harvest of experimen- 
tal observables, with still much more to come, how can the 
GPD information be extracted ? The issue is not trivial as 
we recall that: 



• the Bethe-Heitler (BH) process is another mecha- 
nism which leads to the same final state ep — > ep7 
as DVCS. In the BH process, the final state pho- 
ton is radiated by the incoming or scattered elec- 
tron and not by the nucleon itself. Therefore, the 
BH process, which dominates the cross sections in 
some kinematic regions, carries no information about 
GPDs. However, it is relatively precisely calculable 
in Quantum Electro-Dynamics (QED) given the nu- 
cleon form factors, which are quite precisely known 
at the kinematics we are presently interested in, i.e. 
small t. 

• in the QCD leading twist and leading order approx- 
imation, which is the frame of this study, there are, 
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for DVCS, four independent GPDs: H, E, H and E 
which correspond to the various spin and helicity ori- 
entations of the quark and nucleon in the handbag di- 
agram. These four GPDs depend on three variables 
x, £ and t. Decomposing the DVCS amplitude into 
real and imaginary parts leads to eight GPD-related 
quantities. We will call them the Compton Form Fac- 
tors (CFFs) and they are the quantities which can in 
principle be extracted from DVCS experiments. Fol- 
lowing our conventions introduced in Refs. IU6LU7I1 . 
these eight CFFs are: 

h Re = p [ dx [h(x, e, t) - H(- x , e, t)\ c + (x, ai) 

Jo 

E Re = P [ dx[E(x,C,t) -E(-x,£,t)}C + (x,0(2) 
Jo 

H Re = P f dx \h(x, e, t) + H(-x, £, tj\ C- (x, 0) 
Jo L J 

E Re = pJ dx[E(x,£,t) + E(-x,Z,t)\c-{x,ZU) 

H Im = H(Z,S,t)-H(-S,S,t), (5) 

E Im = E(Z,S,t)-E(-£,S,t), (6) 

Hi m = H(S,S,t)+H(-Z,Z,t) and (7) 

E Im = E(Z,Z,t) + E(-Z,t,t) (8) 



with 



(9) 



In the QCD leading twist and leading order approxima- 
tion, these eight CFFs depend only on £ (or equivalently 
x B ) and t. 

In Refs. iflil [nil , we have developped and applied a 
largely model independent fitting procedure which, at a 
given experimental (x B , —t) kinematic point, takes the 
CFFs as free parameters and extracts them from DVCS ob- 
servables using the well established QCD leading twist and 
leading order DVCS+BH theoretical amplitude. The ex- 
pression of this amplitude can be found, for instance, in 
Ref. OH. With this procedure, we have fitted in Ref. iflill 
the JLab Hall A proton DVCS beam-polarized and unpo- 
larized cross sections. We could then extract the H Irn and 
H Re CFFs at < xb >= 0.36 and for several t values 
with average uncertainties of the order of 35% for Hj m 
and larger for H Re . In Ref. lfl7ll . we have fitted several 
HERMES beam-charge, beam-spin and transversely po- 
larized target spin asymmetries. We could then extract at 
< x R >= 0.09 and for several t values, the same CFFs 
H Irn and H Re with roughly similar uncertainties as for 
JLab. 

The sources of uncertainty in our approach stem, on the 
one hand, from the experimental errors of the data that we 
fit, and on the other hand, from the fact that we take in 



our fits practically all CFFs as free parameters, with rela- 
tively large and conservative bounds. There are therefore 
minimum conjectures and surmises in our work, which is 
certainly highly valuable. However, given that we gener- 
ally fit limited sets of data and observables, our problem 
is in principle underconstrained. The consequence is that 
there are maximum correlations and interferences between 
our fitted parameters, hence the relatively important error 
bars in our results. In these extremely conservative con- 
ditions, it is nevertheless remarkable that we managed, in 
our previous works, to extract several CFFs, at different 
energies, with well-defined uncertainties, fitting the very 
limited available data. The reason for this convergence 
of a few CFFs, in spite of the underconstrained nature of 
the problem, is that some observables are in general dom- 
inated by some particular CFFs, like BSAs by H Im and 
beam charge asymmetries and cross-sections by H Re . Our 
uncertainties can only decrease in the future as, on the one 
hand, larger (and more precise) sets of data and observables 
sensitive to different CFFs become available and, on the 
other hand, theoretical constraints which allow to reduce, 
in the most model-independent way possible, the range of 
variation, or even the number, of the CFFs come forth (for 
instance, dispersion relations Il9l422ll "). 

We mention that related DVCS fitting studies have been 
published this past year j23ll24ll . They resulted also in the 
extractions of the H Im and H Re CFFs with central val- 
ues consistent with ours, although with smaller uncertain- 
ties. These fits are however model dependent. They ei- 
ther neglected all CFFs other than H Im and H Re or (and) 
assumed a functional shape for the CFFs, allowing to fit 
several (i B , —t) points simultaneously. The uncertainty 
associated to the model dependence and hypothesis enter- 
ing these approaches is then very difficult to estimate and 
to take into account properly. Nevertheless, each of these 
model-dependent and -independent approaches have their 
own merits and values. The fact that they all result in con- 
sistent and compatible central values for the fitted H Im 
and H Re CFFs gives mutual support and credit for each 
of them. 

In this letter, we continue our model-independent fitting 
approach focusing this time on the CLAS BSAs and ITSAs, 
which we had not considered so far. As was already men- 
tioned in our pioneer work on the subject fltHl . fitting only 
the CLAS BSAs, without any model-dependent hypothe- 
sis or input, was not constraining enough. In other words, 
fitting only BSAs, i.e. only one observable, with seven 
or eight unconstrained CFFs does not lead to well-defined 
solutions. However, inspired by our recent rather succes- 
ful fit of a series of HERMES asymmetries 11711 . we now 
want to investigate if, with the addition of a new observ- 
able, namely the ITSAs, to be fitted simultaneously with 
the BSAs, progress can be made. Indeed, ITSAs have been 
measured by the CLAS collaboration as well and have ac- 
tually received little attention from GPD phenomenologists 
so far. The BSAs being in general dominantly sensitive to 
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Him an d the ITS As to Hi m IU6LI25I1 . our expectation is to 
extract some quantitative information on these two particu- 
lar CFFs, which would be brand new information for Hi m 
in particular. 

Let us describe these CLAS data. Regarding BSAs, the 
Hall B collaboration has measured their (f) distribution at 
57 (x B , —t, Q 2 ) points {(f) is the standard angle between 
the leptonic and hadronic planes of the DVCS process). 
The values of the x B variable extend from ~0.13 up to 
~0.46, those of the —t variable from ~0.13 GeV 2 up to 
wl.3 GeV 2 and those of the Q 2 variable from GeV 2 
up to ~3.3 GeV 2 . The amplitude of these BSAs range from 
«0upto«0.3 m. 

Regarding ITSAs, the data is much more scarce: in 
Ref. II 1211 . only their sin(</>) moment is available at a few 
(x B , —t, Q 2 ) points. The average kinematics of this 
whole set of data is < £ >=0.16 (i.e. < x B >«0.275), 
< -t >=0.31 GeV 2 and < Q 2 >=1.82 GeV 2 . Within 
this phase space, the 1TSA sin(</>) moments (which we 
will designate as A^f) following, for instance, Ref. H3l) 
have been extracted differentially, either for three x B val- 
ues (psO.20, 0.29 and 0.40) or for three -t values («0.15, 

0. 24 and 0.43 GeV 2 ). There are therefore six ITSAs avail- 
able, which are actually not statistically independent since 
they are issued from the same set of data which has been 
binned either in x B or in t. However, even if scarce, these 
data are extremely valuable as we will show in the follow- 
ing. Their amplitudes range from ~ 0.07 to ~ 0.38, with 
uncertainties extending from 30% to more than 100%. 

Following our notation of Ref. [17], the BSAs can also 
be denoted as A LU and the ITSAs moments as A s ^j . We 
will use this notation on our figures. 

In the first stage, we will see what information one can 
extract from the simultaneous fit of one 1TSA and one BSA 
whose (x B , —t, Q 2 ) values approximatively match. In the 
second stage, we will fit simultaneously one 1TSA and the 
two or three BSAs which have the same (x B , —t) values, ir- 
respective of the Q 2 value. Assuming the dominance of the 
QCD leading twist and leading order DVCS contribution, 

1. e. that CFFs do not depend on Q 2 , the idea is to improve 
the statistical accuracy and increase the constraints on our 
fitting procedure. 

Among the six available ITSAs, we first focus on the 
ITS A which has been measured at the kinematic point (x B , 
-t, Q 2 ) =(0.29,0.31,1.82). Unfortunately, the BSAs and 
the ITSAs are issued from two different Hall B experi- 
ments so that they have not been measured at exactly the 
same average kinematics. The matching of the kinematics 
between the different observables can thus be only approx- 
imate. Among the 57 BSAs, the two BSAs whose kinemat- 
ics is the closest of the 1TSA kinematics that we focus on, 
are at (0.25,0.28,1.69) and (0.25,0.28,1.95). We therefore 
note at this stage the differences between the x B values 
(0.25 vs 0.29), the -t values (0.28 vs 0.31) and Q 2 val- 
ues (1.69 or 1.95 vs 1.82). There is a third BSA which 



has the same x B and t values as the two BSAs just men- 
tioned but whose Q 2 is equal to 2.21 GeV 2 . We will use 
this extra BSA in our Q 2 -independent "second stage" fit- 
ting. We recapitulate these four kinematic points on which 
we presently focus in Table U 





< x B > 


< -t > 


<Q Z > 


1TSA 


0.29 


0.31 


1.82 


BSA1 


0.25 


0.28 


1.69 


BSA2 


0.25 


0.28 


1.95 


BSA3 


0.25 


0.28 


2.21 



TABLE I: Summary of the four kinematical points which have 
approximately the same {x B , —t) values. We will fit four 
"topologies": (1TSA+BSA1), (1TSA+BSA2), (1TSA+BSA3) and 
(1TSA+BSA1+BSA2+BSA3). 



For the fitting procedure, as in Refs [16, 17], we mini- 
mize our theoretical calculations of the DVCS observables 
based on the well-known QCD leading twist and leading 
order DVCS+BH amplitude, by the standard x 2 function, 
using MENUIT [26]. We recall that the parameters to be fit- 
ted are the CFFs of EqsQJSJ As in Refs Illq,ll7ll . we have 
actually considered only seven CFFs, setting E Im to zero. 
This is based on the theoretical guidance which approxi- 
mates the E GPD by the pion exchange in the t-channel 
whose amplitude is real. With the hypothesis of the dom- 
inance of the leading twist amplitude of the DVCS pro- 
cess, this is the only model-dependent assumption that en- 
ters our fitting procedure. A last feature entering our fitting 
process is that we have to bound the domain of variation 
of the fitting parameters. Without bounds, our fits which 
are in general underconstrained, would probably not con- 
verge and/or would yield values for the fitted parameters 
with infinite uncertainties. Following what we have done 
and explained in details in Refs. fl6L \Yfa . we bound the 



allowed range of variation of the CFFs to ±5 times some 
"reference" VGG CFFs. VGG 01 H3] is a well-known 
and widely used model which provides an acceptable first 
approximation of the CFFs, as shown in our previous stud- 
ies 

010 

and as will be confirmed furtherdown in the 
present work. We do not really consider this as a model- 
dependent input since this allowed deviation of a factor ±5 
with respect to the VGG model values is extremely conser- 
vative. We recall that GPDs have to satisfy a certain num- 
ber of normalization constraints in general, these being all 
fulfilled by the VGG model. Finally, the problem at stake 
being non-linear and the parameters being correlated, we 
use MINOS for the uncertainty calculation on the resulting 
fitted parameters Il26ll . 

Before presenting our results, we also want to outline the 
point that our aim is to fit the 1TSA and BSA(s) with the 
same CFFs, which therefore should correspond to unique 
(x B , —t) values. As the data of Table U do not, unfortu- 
nately, have exactly the same kinematics, as mentioned pre- 
viously, there is an ambiguity in defining the precise (x B , 
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— t) values of the fitted CFFs. We will consider that the 
CFFs that we will fit to the kinematic points of Table U cor- 
respond to the values of the BSA kinematics, i.e. (0.25, 
0.28), as the BSA observable is in general the most signif- 
icant statistically. There is clearly an approximation here, 
which we will make for the moment, lacking better solu- 
tion, in order to make progress. This approximation is to 
some extent supported by the VGG model which predicts 
about 8% difference for H Im between Xg=0.25 and 0.29 
(at — 1=0.28 GeV 2 ) and less than 2% for H Im for this same 
kinematic change. We stress that this kinematical matching 
problem is of a rather trivial nature and it is sufficient that 
future experiments, measuring different observables, sim- 
ply agree to analyze data at the same central kinematics to 
avoid this extrapolation issue. 

We now display in Fig. [2] the result of our fits. The 
left panel shows the eft distribution of the three BSAs men- 
tionned above and the right panel the sin <fi moment of the 
1TSA. The dashed curves are the results of the fit of the 
1TSA (of the right panel) with each individual BSA. The 
thick solid curves are the result of the fit of this same 1TSA 
with the three BSAs simultaneously (these three BSAs hav- 
ing the same —t) values but different Q 2 values, see 
TableU). On the right panel, the four empty circles show the 
corresponding results of the fit for the 1TSA: the first three 
for the fit with the individual BSAs and the fourth one for 
the fit with the three BSAs simultaneously. For compari- 
son, we also show in this figure the predicted results for the 
BSAs and the 1TSA of the standard VGG model HHStI]. It 
is seen that the VGG model overestimates the three BSAs 
by approximately 0.1 (i.e. ~ 30%) and underestimates the 
1TSA by roughly the same proportion. 

We now show in Fig.|3]the fitted values, with their error 
bars, of the only two CFFs, H Im and H Im , out of seven, 
that came out of our fitting procedure with finite MINOS 
uncertainties. We recall that the MINOS uncertainties cor- 
respond to a deviation of 1 from the value which minimizes 
\ 2 . These uncertainties can be asymmetric if the x 2 func- 
tion is not symmetric around the minimum, which is the 
sign of a non-linear problem in general. The fact that only 
H Im and H Im converge in our fitting process reflects, as 
was mentioned earlier, the particular sensitivity of the BSA 
and 1TSA observables, respectively, to these two CFFs. 
The other five CFFs did not converge in our fitting pro- 
cedure to some well defined value or domain: either their 
central value reached the boundaries of the allowed domain 
of variation or MINOS could not reach the x 2+i value to 
fully determine the associated uncertainties. These features 
were well studied lfl6l[l7ll in our earlier works. They reflect 
the fact that the contribution to the \ 2 °f those CFFs which 
didnotconverge is relatively weak and that the fit is barely 
sensitive to them. However, it is important to include them 
in our fit because they play a role, through correlations, in 
the determination of the error bars on the two "convergent" 
CFFs. 
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FIG. 2: Comparison of our fit results with the experimental data. 
The three left panels show the three experimental BSAs (i.e. 
Ajju) as measured in Ref. 11(31 (solid circles). The right panel 
shows the experimental 1TSA moment (i.e. A s ^y)as measured 
in Ref. [ 12] (solid circle). The four panels correspond to the four 
(xg, —t, Q 2 ) kinematic points presented in Table U (from left to 
right: BSA1, BSA2, BSA3 and 1TSA). All four observables have 
approximately the same (xg, —t) values, taken as (0.25,0.28), 
but different Q 2 values. On the BSA panels, the dashed solid 
curve is the result of our fit, fitting only the individual BSA of 
the relevant panel with the 1TSA, i.e. from left to right, fit of 
(BSA1+1TSA), (BSA2+1TSA) and (BSA3+1TSA). The solid line 
is the result of our fit, fitting simultaneously the three BSAs and 
the 1TSA, i.e. (B S A 1 +B S A2+B S A3+1TS A) . This latter fit there- 
fore assumes that CFFs do not depend on Q 2 . The fit results of 
these four "topologies" for the 1TSA are displayed in the right 
panel. On the BSA panels, the dotted curve is the prediction of 
the standard VGG model. Its prediction for the 1TSA is displayed 
as the empty cross. 



In Fig. [3] we display four sets of results for H Im and 
H Irn , which correspond to the four "topologies" that we 
mentionned earlier: three sets correspond to the fits of the 
1TSA with each one of the BSAs at Q 2 =1.69, 1.95 and 2.21 
GeV 2 and the fourth set (in the box in Fig. [3]) corresponds 
to the simultaneous fit of the 1TSA with the three BSAs. In 
this latter case, the underlying assumption is that CFFs do 
not depend on Q 2 . As could be expected, the resulting un- 
certainties are smaller for the both CFFs in this latter con- 
figuration, as more statistics and contraints enter into play. 
We observe that all four configurations yield compatible 
results within error bars, which are between 25% and 50% 
in average. The simultaneous fit of the three BSAs and of 
the 1TSA yields an approximate average of the fits using 
only one BSA and the 1TSA. From the uncertainties on the 
CFFs that we obtain, it is clear that no QCD evolution or 
twist effect can be discerned. It then seems reasonnable to 
fit simultaneously observables at the (approximately) same 
(xg, —t) points and different Q 2 values. 

In Fig. [3l we have also displayed, for each of the four 
fit topologies, two results, aimed at illustrating the depen- 
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FIG. 3: Results of our fits to the data displayed in Fig. [3] 
for the Hj m and Hi m CFFs. The empty squares (circles) 
show our results when the boundary values of the domain 
over which the CFFs are allowed to vary is 5 (3) times the 
VGG reference values. Both results have been slightly shifted 
left (right) from the central Q 2 value for sake of visibil- 
ity. The four sets of results correspond, from left to right, 
to the fits (1TSA+BSA1), (1TSA+BSA2), (1TSA+BSA3) and 
(ITS A+BS A 1 +B S A2+B S A3) as indicated in Table U In par- 
ticular, the (1TSA+BSA1+BSA2+BSA3) fit, which corresponds 
to an average Q 2 value of 1.95 GeV 2 , is displayed within a 
box. The empty diamond indicates the results of our fits, in the 
(1TSA+BSA1+BSA2+BSA3) "topology", when only the H and 
H GPDs are taken as fitting parameters. The empty cross indi- 
cates the VGG prediction. 



dence of our results on the boundaries of the domain of 
variation allowed for the CFFs. The empty squares show 
our results when the CFFs are limited to vary within ±5 
times the VGG reference values while the empty circles 
shows these results for boundaries equal to ±3 times these 
same VGG reference values. Of course, the smaller the do- 
main of variation, the smaller the uncertainties on the fitted 
CFFs. This shows the overall stability and robustness of 
our fitting process since the values of these boundaries do 
not affect strongly the central values of the fitted CFFs. We 
also checked that the fit results were not dependent on the 
precise starting values of the CFFs when we begin our fit: 
irrespective of the starting values, the minimization would 
essentially always converge to the same central values and 
uncertainties for H Im and H Im . 

We further show in Fig.[3]the result of our fit if we take as 
fitting parameters only the H and H GPDs, i.e. four CFFs 
(H Re , H Re , H Im and H Im ), instead of seven. For this con- 
figuration, we have fitted the three BSAs and the 1TSA si- 
multaneously. The central values for H Im and H Im are in 
very good agreement with the ones previously determined 
when all CFFs were taken into account (with, though, a 
some slight decrease of the central value of H Im ). The 
obvious difference is that, as could be expected, the associ- 



ated uncertainties are smaller, particularly for H Im . There 
is not too much effect for H Im . This can probably be at- 
tributed to the fact that when only two GPDs enter the fit, 
the main source of uncertainty comes from the statistics 
of the observables to be fitted and no more from the cor- 
relations between the fitting parameters. Indeed, H Im is 
mostly sensitive to the three BSAs (which are simultane- 
ously fitted) with each having smaller errors than the 1TSA, 
while H Im is mostly sensitive to the 1TSA which has a ~ 
25% error bar. We do not display the comparison of these 
"2 CFFs" fit with the data in Fig.|2]in order not to overload 
the figure, but the x 2 1S equally good to the fits with all 
CFFs. In the latter case, the normalized x 2 is found to be 
equal to 1.27 while in the former case the normalized x 2 
is 1.16. These good results obtained when fitting with only 
H Im and H Im mean that it is indeed possible to correctly 
fit the data with only these two GPDs instead of four. This 
however does not mean that this is the true solution and 
that the other GPDs should consequently be ignored or ne- 
glected. The large error bars that we obtain when all GPDs 
precisely reflect this lack of knowledge on the other CFFs: 
our uncertainties incorporate all our ignorance about the 
other GPDs and all their full potential influence. 

We finally display in Fig. [3] the predicted values of the 
corresponding VGG CFFs (empty crosses), which are Q 2 
independent. It is noted that the VGG H Im tends to lie 
above the fitted H Im while the VGG H Im tends to lie be- 
low the fitted H Im . This is a straightforward reflection of 
what was observed in Fig. [2] where the VGG BSAs curves 
were overestimating the data while the VGG 1TSA point 
was underestimating the data. The overestimation of the 
VGG H Im , with respect to the fitted central value, was also 
observed in our study of the HERMES data I17I1 . 

We have so far focused on the particular 1TSA measured 
at x B =0.29 in order to establish and understand the ba- 
sic features and results of our approach. We now turn to 
the 1TSA measured at another x B value, i.e. x B =0.40, 
for which there are several BSAs which have neighboring 
(xb, —t) values. These data points are indicated in Ta- 
ble |IIJ This time, none of the Q 2 values match each other 
and there is also a more significant difference between the 
Xb values. Supported by our previous study which showed 
that the simultaneous fits of several observables at differ- 
ent Q 2 's appeared to converge to some average of individ- 
ual Q 2 fits, we attempt to fit simultaneously the four data 
points (i.e. 3 BSAs and 1 1TSA) of Table M We men- 
tion the VGG predictions: 13% change for H Im between 
x B =034 and x B =0.40 (for -t=0.30 GeV 2 ) and 7% change 
for H Im for the same kinematical variation. In front of the 
anticipated ~30% error bars to be issued from our fits, it is 
not unreasonnable to neglect, in a first approach, this small 
Xb variation. 

Within this approximation, we are again able to ex- 
tract values for the two CFFs H Im and H Im which we 
then consider to correspond to the BSAs' kinematics (xb, 
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1TSA 


0.40 


0.31 


1.82 


BSA1 


0.34 


0.30 


2.3 


BSA2 


0.34 


0.28 


2.63 


BSA3 


0.35 


0.29 


2.97 



TABLE II: Selection of kinematic points measured by the CLAS 
collaboration which have approximatively the same (ig, — t) val- 
ues, around a;s=0.35. 

— i)=(0.35, 0.29), based on the statistical dominance of the 
BSAs. We show the resulting values of H Im and H Im in 
Table HTI1 along with the values we obtained previously for 
these two CFFs at x B =0.25 when we fitted simultaneously 
all points of Table U (i.e. values of the data points in the 
"box" of Fig. |3). We observe, although error bars are not 
negligible, the general trend that, at fixed t (~0.28 GeV 2 ), 
H Im tends to increase, as Xb goes from 0.35 to 0.25, while 
H Im remains rather constant. 



< x B > 


Him 


Him 


0.25 


i.oo_ 71 




0.35 


U - DZ -0.18 


63+"' bu 



TABLE III: Results of our fits for the Hi m and Hi m CFFs from 
the CLAS BSAs and 1TSA, at fixed t Ri0.28 GeV 2 , for two dif- 
ferent x b values. 

We recall that we were able in earlier work to also ex- 
tract values for Hi m at different x B values and at almost 
the t values considered here (~-0.28 GeV 2 ). For recall, 
in Ref. |fl^l . we fitted the JLab Hall A data which have 

< x B >«0.36 and in Ref 03], we fitted the HERMES 
data which have < x B >~ 0.09. While the JLab Hall 
A data were taken precisely at < —t >=0.28 GeV 2 , the 
HERMES data were given for < -t >=0.20 GeV 2 and 

< —t >=0.42 GeV 2 . In a very simplistic way, we decide 
to interpolate between these two —t values by simply av- 
eraging our fitted H Jm CFFs at these two —t values. We 
also average quadratically the positive and negative error 
bars. We thus end up with some average HERMES Hi m 
CFF at < -t >~ 0.30 GeV 2 and < x B >~ 0.09. We 
can then obtain a x B dependence of our fitted Hi m 's using 
our JLab and HERMES analysis results. Fig. [4] shows this 
x B dependence, compiling our results from the indepen- 
dent analysis of the JLab Hall A data, the HERMES data 
and the presently analyzed CLAS data. 

It turns out that around Xb=035, both JLab Hall A and 
CLAS data are available. It is comforting to note the de- 
cent agreement, within error bars, of the two extracted val- 
ues for Hi m . True, the error bars are not small and it might 
appear not so challenging to have an agreement with such 
uncertainties. Nevertheless, we would rather take this as- 
pect as support for our realistic evaluation of the error bars 
on our fitted CFFs. 

In order to illustrate this point, we also plot in Fig. [4] 



(open diamond around x B =0.35) our fit results of the 
CLAS BSAs and 1TSA at x B =035 with only H and H 
as fitting parameters, i.e. four CFFs instead of seven. It 
is seen that the central value result which, as could be ex- 
pected, has a quite smaller error bar, is slightly shifted with 
respect to the central value result when the seven CFFs are 
taken into account. However, both results remains well 
compatible within error bars. We did the same exercice 
with the JLab Hall A data at x B =0.36. In Fig. 01 on the 
one hand, the open triangle shows the result of the fit to 
the JLab Hall A unpolarized and beam-polarized cross sec- 
tions at < —t 0.28 GeV 2 , using only the H GPD, i.e. 
two CFFs. On the other hand, the open diamond around 
Xs=0.36 the same fit using only H and H, i.e. four CFFs. 
Let us mention that these two fits give x 2 values very close 
to one and describe very well the data points. Since these 
two fits use only a few CFFs, the error bars on the resulting 
H Im 's are quite (artificially ?) small compared to the error 
bar resulting from the fit with seven CFFs (open square). 
However, it is striking to see how the central value of H Im 
can shift. With only H entering the JLab Hall A data fit, 
Him (open triangle near x B =0.36) is around two, at the 
top of the error bar of the seven CFFs fit result (open square 
near x #=0.36). This slight increase of Hi m is actually very 
consistent with the result of Ref. Il23ll which used the same 
assumption of H dominance. However, this result is now 
clearly inconsistent with the result for Hi m issued from our 
fit of the CLAS data (open square near x B =035). 

Turning to the configuration where both H and H enter 
the fit of the same JLab Hall A data, Hi m has now dramat- 
ically dropped by a factor more than two (open diamond 
near x B =0.36). However, it has now actually become com- 
patible with the result issued from the CLAS data, either 
open square or diamond near Xb=035. Indeed, note that 
both CLAS Him results, i.e. issued from the fit with only 
H on the one hand and with the seven CFFs on the other 
hand, are compatible within error bar. This is a clear-cut 
illustration of the meaning of the large error bars which 
result from our fits when we use seven CFFs: they cau- 
tiously and realistically reflect the underconstrained nature 
of our problem, i.e. fitting only a couple of observables 
with many parameters (i.e. the CFFs), and all the variation 
and potential weight of these parameters. 

To summarize this discussion, in the framework of our 
analysis (i.e. leading twist and leading order QCD and 
the few kinematic approximations mentionned earlier), it 
doesn't seem to us possible to find a consistent value of 
Him to fit both the JLab Hall A and CLAS data if only H 
enters the fit. The minimum scenario seems that H be in- 
cluded, the ultimate one being of course that all CFFs be 
included. It is interesting to mention that Ref. ll24ll reached 
some similar conclusion in a model-dependent approach, 
confirming the hint that GPDs other than H (and possibly 
H) do play a significant role at the JLab kinematics. Let 
us stress again that at this kinematic point, < Xb >~0.35, 
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our values of H Im were determined by the fitting of inde- 
pendent DVCS experiments, i.e. JLab Hall A and CLAS, 
and rather different observables: polarized and unpolarized 
cross sections for the JLab Hall A analysis and BSA and 
1TSA for the present CLAS analysis. Although beam po- 
larized observables are common to the two experiments, it 
is encouraging to observe that different paths can lead to 
consistent results, as it should be. 

Now, more generally, taking into account the HERMES 
data, we observe in Fig.|4]that the general tendency is that, 
at fixed t, H Im increases with decreasing x B . This is rem- 
iniscent of the x b dependence of standard partem distribu- 
tions. The VGG prediction is also shown in Fig. [4] and, 
although it overestimates most of the fitted central values, 
it displays the same behavior. 



same x B and —t values with, however, different Q 2 val- 
ues. We list those points in table ||Vj Comforted by our 
reasonable results presented in Fig. [3] we adopt the same 
approach and fit, with the seven CFFs as fitting parame- 
ters, simultaneously the three BSAs and the 1TSA at each 
of the three (x B , —t) points of Table [IV] which all have 
a common i B value (i.e. ~ 0.25). Again, only the H Im 
and H Im CFFs systematically come out from our fits with 
finite error bars. Fig.[5]shows our results and reveals the t- 
dependence (at Xb ~0.25) of the H Im and H Irn CFFs. We 
again display for each t value two results corresponding, 
like in Fig. [3] to different boundary values for the domain 
of variation allowed for the CFFs, i.e. ±5 (empty squares) 
and ±3 (empty circles) times the VGG reference values. 
We also show in this figure the VGG predictions (empty 
crosses). 



fit of 
HERMES 



□ 7 CFFs 
only H,H 
A only H 



fit o-' 
JLab Hall A 



t 
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fits of JLab Hall B 

0.1b 0.2 0.2b 0,3 0.35 0.4 0.4b 



FIG. 4: x B dependence at fixed -i=0.28 GeV 2 of the fitted 
Hj m (empty squares) according to our analyzes of the JLab 
Hall A data [3] « x B >=0.36), of the HERMES data 01 
(< xb >=0.09) and of the present analysis (< xb >=0.25 and 
< xb >~0.35), using the seven CFFs. The open diamond, 
slightly shifted left (for visibility) of the open square point at 
x B «0.35 is the result of the CLAS BSAs and ITS A at x B -0.35 
using only the H and H GPDs. The open triangle, slightly shifted 
right (for visibility) of the open square point at xb ~0.36 is the 
result of the JLab Hall A unpolarized and beam-polarized cross 
sections at ££=0.36 using only the H GPD. The open diamond, 
slightly shifted right (for visibility) of the open square point at 
xb ~0.36 is the result of the JLab Hall A unpolarized and beam- 
polarized cross sections at 2:^=0.36 using only the H GPD and 
H GPDs. The empty cross indicates the VGG prediction. 

We now finally turn our attention to the t dependence 
of the ITSAs. The CLAS collaboration has extracted the 
ITSAs at fixed x B (^0.25) for three different < —t > val- 
ues: 0.15, 0.24 and 0.43 GeV 2 . For each of these -tpoints, 
we can identify three BSAs which have approximately the 





<x B > 


< -t > 


<Q A > 


1TSA 


0.27 


0.15 


1.82 


BSA1 


0.24 


0.15 


1.65 


BSA2 


0.24 


0.14 


1.89 


BSA3 


0.25 


0.14 


2.16 


1TSA 


0.27 


0.24 


1.82 


BSA1 


0.25 


0.28 


1.69 


BSA2 


0.25 


0.28 


1.95 


BSA3 


0.25 


0.28 


2.21 


1TSA 


0.27 


0.43 


1.82 


BSA1 


0.25 


0.49 


1.70 


BSA2 


0.25 


0.49 


1.95 


BSA3 


0.25 


0.49 


2.20 



TABLE IV: Selection, for our i-dependence study, of the three 
(xb, —t) kinematic points measured by the CLAS collaboration 
which have 1 1TSA and 3 BSAs at approximately the same ig 
values (« 0.25). 

In Fig.|5j regarding H Im , we note a smooth and typical 
fall-off with —t which was also observed in our previous 
JLab Hall A and HERMES studies |E El- The figure 
also confirms that the standard VGG parametrisation, in 
general, overestimates our fitted values. This is particu- 
larly the case at low t (this was also observed at HERMES 
energies U3l). Regarding H Im , although the uncertainties 
are large, the /j-slope appears to be much less pronounced 
and it even seems that there is a drop towards as t goes 
to (although a constant and flat i-dependence can also be 
in order within error bars). We find again that VGG under- 
estimates this CFF, in particular as | t | grows. Overall, 
the VGG /j-slope is markedly different from the one of the 
fitted Hi m . 

Finally, the diamonds in Fig. [5] show the results of our 
fit when only the H and H GPDs are taken as fitting pa- 
rameters, i.e. setting to all other GPDs. We observe the 
same features as previously (see Fig. 0. The central val- 
ues are in very good agreement with the ones determined 



8 







H Im (x B =.25,t) 


H Im (x B =.25,t) 


: * 
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0.2 0.4 0.2 0.4 0.6 



-t (GeV J ) 



FIG. 5: ^dependence of our fitted Hj m and CFFs at the 
kinematic points of Table |IV| The empty squares (circles) show 
our results when the boundary values of the domain over which 
the CFFs are allowed to vary is 5 (3) times the VGG reference 
values. The empty diamonds indicate the results of our fits when 
only the H and H GPDs are taken as fitting parameters. The 
empty crosses indicate the VGG prediction. At the lowest t value 
of Hj m , the empty cross and diamond happen to be superim- 
posed. 

from the fits in which all CFFs were included (though with 
a systematic decrease by ~ of 15% of the central value of 
H Irn ). The main effect is to reduce the uncertainties on the 
fitted CFFs: very strongly for H Im and only slightly for 

Hl m . 

To summarize this work, we have analyzed the beam 
spin asymmetry and the longitudinally polarized target spin 
asymmetry of the Deep Virtual Compton Scattering pro- 
cess recently measured by the CLAS collaboration. We 
have used a fitter code, largely model-independent, based 
on the QCD leading-twist and leading order DVCS+BH 
amplitude, which takes as fitting parameters GPD CFFs. 
Even though we fit only two asymmetry observables with 
seven CFFs, two CFFs, H Im and H Im , come out systemat- 
ically from our fits with stable and well defined central val- 
ues and uncertainties (of the order of 30% in average). The 
reason is that the two observables we fit are well known to 
be dominantly sensitive to these two CFFs. It is worth not- 
ing that with only BSAs to fit, there is no convergence of 
our fits, while with the addition of a single observable, i.e. 
the 1TSA, solutions become relatively well defined. 

In this work, a few approximations have been made, 
mostly due to the present lack of sufficiently precise and 
numerous data. We recall that only six ITSAs were avail- 
able to us: for three —t values at fixed x B and Q 2 and for 
three Xb values at fixed —t and Q 2 values. The approx- 
imations that we did were to simultaneously fit BSAs and 
ITSAs taken at slightly different x B and —t values and also 
at different Q 2 values. This latter approximation is, in any 
case, along the line of the main starting assumption of this 



work: the dominance of the QCD leading twist and lead- 
ing order of the DVCS amplitude. In those conditions, we 
have been able to determine the s B - and t- dependences of 
the H Irn and H Im CFFs (respectively at fixed t and fixed 
Xb)- In particular, we put in evidence a much flatter t- 
dependence for H Im than for H Im . We also illustrated, 
by comparing our fits at roughly the same kinematics of 
the independent JLab Hall A and CLAS data, the impor- 
tance and, even the necessity, of taking into account several 
GPDs in order to obtain compatible results. 

While there have lately been a couple of other works 
aiming at fitting DVCS data and extracting H Im , this is the 
first one allowing access to H Im and, in a largely model- 
independent way, determine some first numerical value for 
it. The "price" to pay for our model-independency is that 
we obtain relatively large uncertainties. Several DVCS ex- 
periments aiming at measuring more precisely the observ- 
ables analyzed in this work and also aimed at measuring 
new observables, such as transversely polarized target spin 
asymmetries and cross sections are under way in the near 
future. We expect our fitting techniques to be more and 
more fruitful and efficient as these new precise and nu- 
merous data become available, along with theoretical GPD 
modelling progress which can reduce the domain of varia- 
tion of the fitted CFFs or their number. 
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